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Abstract 



The second-quantization of a scalar field in an open cavity is formulated, 
from first principles, in terms of the quasinormal modes (QNMs), which are 
the eigensolutions of the evolution equation that decay exponentially in time 
as energy leaks to the outside. This formulation provides a description involv- 
ing the cavity degrees of freedom only, with the outside acting as a (thermal 
or driven) source. Thermal correlation functions and cavity Feynman propa- 
gators are thus expressed in terms of the QNMs, labeled by a discrete index 
rather than a continuous momentum. Single-resonance domination of the 
density of states and the spontaneous decay rate is then given a proper foun- 
dation. This is a first essential step towards the application of QNMs to cavity 
QED phenomena, to be reported elsewhere. 
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I. INTRODUCTION 



In this paper we are concerned with quantum fields in open cavities — the obvious 
example of ultimate interest would be electromagnetic fields in optical cavities and the 
associated problem of cavity quantum electrodynamics (CQED). Such systems are open 
because energy leaks to the outside (e.g., via output coupling), and as dissipative systems 
cannot normally be quantized on their own rather, one must also consider the bath into 
which energy escapes, so that the total universe is conservative. Thus such cavities, say of 
linear dimension a, can be embedded in a universe of dimension A ^ oo. One can then 
quantize on the modes of the universe, which are labeled by a wavenumber p spaced by 
Ap ~ vr/A ^ 0. The field quanta are then created or destroyed by operators a^(p) and a{p), 
and higher-order processes involve integrals Jdp ■ ■ ■ . 

On the other hand, these cavities often have a very small amount of leakage, characterized 
by a parameter e = -C 1, where the quality factor of the cavity can be as high as Q ~ 10^ 
or more. If this is the case, the intuition developed from a closed cavity, also of length a, 
should be relevant. A closed cavity is a conservative system, with normal modes (NMs) 
labeled by a discrete index j = 1,2, ■ ■ -, where the wavenumber is pj ~ jvr/a, Ap ~ ir/a. 
Field quanta in such a closed cavity are created and destroyed by operators a] and aj, and 
higher-order processes involve discrete sums J2j ■ ■ ■• Can quantum fields in an open cavity be 
described in a similar way — in terms of discrete modes and the corresponding operators? 
If this is possible, computations will be simplified and will correspond to physical intuition, 
with each term j associated with a cavity "mode" . The connection with the limit of a closed 
cavity (e —>■ 0) would also become manifest. 

Quantization of a closed system relies on its NMs; the counterparts in an open system 
are the quasinormal modes (QNMs), which are again factorized solutions 

<f>{x,t) = fj{x)e-'^^\ (1.1) 

where Imuj < because of the loss of energy. Each QNM corresponds to a resonance, 
with a width 7^ = | Imc(jj|. The purpose of this paper, in short, is to develop a formalism 
whereby field quantization can be implemented in terms of these QNMs, and to define and 
study operators a], aj for these modes. Specifically, one wishes to express field correlation 
functions, Feynman propagators and other quantities in terms of the QNMs. The dissipative 
nature of the system is then contained in the QNMs themselves. 

The advantages for CQED would be obvious. The simplest phenomenon to which such a 
formalism can apply is the well-known enhancement (or suppression) of spontaneous decay 
rates when the emitted radiation falls on (or between) resonances f^. Resonance domination 
of these processes has been discussed starting with the heuristic argument due to Purcell . 
He proposed that the Fermi golden rule should be generalized: the density of states per unit 
volume, instead of the usual do^uj) = u'^ /{tc'^c^) for vacuum (where c is the velocity of light), 
is to be replaced by d^uj) ~ D/{2'yV) for a D-io\d degenerate QNM of width 7 in a cavity 
of volume V. This leads to an enhancement factor of K = d/do {1/8tt)DQ{X^/V) for 
spontaneous emission on resonance, where A is the wavelength of light emitted and Q is the 
quality factor of the cavity. The essence of this argument is that each resonance counts as 
one state, i.e., in a suitable sense it carries unit weight. While intuitively plausible, this 
statement is difficult to justify formally — since the entire concept of a state, i.e., an NM, 
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falls apart in an open system. However, this argument, and its many variants and extensions, 
would find natural expression in a framework that quantizes on the QNMs, and we shall in 
particular show explicitly below that each resonance carries unit weight. 

In Section |I|, the QNM expansion of classical fields outgoing from a cavity is reviewed. 
The classical results may be organized into two levels. First, under conditions to be specified, 
the Green's function G can be expanded in terms of QNMs. Second, one can try to expand 
the outgoing classical fields (p in terms of QNMs, and to establish a linear space structure 
similar to that for conservative systems. In order to do so, it turns out to be necessary to 
make use of a two-component formalism, dealing with (j) and the conjugate momentum 
at the same time. The linear space formalism is more elegant, but in its simplest form is 
limited to only 1 d. 

The quantum formalism can likewise be approached in two ways. In the first, which we 
shall call the Green's function method (Section |ITTD , one focuses on c-number correlation 
functions and propagators without explicitly expanding the field operator in terms of 
QNMs. The key idea is that the retarded propagator of the quantum theory is exactly 
the same as the classical Green's function G, and the latter has a QNM expansion. Once 
G^ is obtained, it is straightforward to derive a QNM expansion for the correlation function 
F as well. 

In the second, which we shall call the field expansion method, one tries to establish an 
expansion of the quantum field (p, in parallel with the linear space structure established for 
classical fields. This allows us to interpret the expansion coefficients a] and aj, roughly 
speaking, as generalized creation and annihilation operators for the discrete QNMs j. How- 
ever, before doing so, it has to be recognized that quantum fields cannot be constrained by 
the outgoing wave condition — for the simple reason that zero-point (and thermal) fiuctua- 
tions must contain an incoming component. Thus, the first step in developing this method, 
presented in Section is to generalize the field expansion to handle incoming waves as 
well. With this generalization, one can then subject the fields to canonical quantization 
in Section 0. This is done by starting with the universe, a closed hermitian system for 
which the quantization is unambiguously defined. Then, in parallel with the usual removal 
of bath oscillators 0, the outside degrees of freedom are eliminated from the equations 
of motion The results will be equations of motion and commutation relations for the 
discrete operators a] and aj, in which the effects of the outside bath are clearly displayed: 
the loss of energy of each mode by leakage, and the pumping of each mode by the thermal 
or quantum fiuctuations from the outside. 

The formalism is then used to evaluate the correlation function F in Section 0, and the 
results are compared with those obtained from the Green's function method. Interestingly, 
the results appear to be different — those derived from the Green's function method contain 
a single sum J2j over the QNMs, while the field expansion method yields a double sum X^jfc 
with off-diagonal terms. The two are, however, shown to be equal through an identity on G^. 
Recalling that the expansion of the classical field is unique only when the second component (j) 
is considered at the same time, we next show that the expansion of the correlation function 
F is also unique if we consider as well, giving the non-diagonal form. The density of 
states d, which is intimately related to the correlation function, is also expressed in terms 
of QNMs; in particular it is shown that up to corrections of 0{Q~^), each resonance carries 
unit weight in the density of states. 
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The results on the correlation function are then used, in Section [VI1| , to evaluate and 
discuss the Feynman propagator G^, which is the fundamental building block for CQED. 
Again, equivalent diagonal and non-diagonal forms are obtained. Particular attention is 
paid to the equal space propagator G^{x,x,uj), whose imaginary part is related to the life- 
time of an excited atom placed at x. This quantity is discussed in the approximation of 
domination by a single resonance, providing justification for Purcell's heuristic argument 0] 
on the enhancement of spontaneous decay rates. The advantage of using the non-diagonal 
expression is again emphasized. 

In Section |V11J| , a very simple example is studied explicitly, and its correlation function 
and energy density are expressed in terms of a sum over QNM contributions. 

Some final remarks are then given in Section |X|. We stress that this paper is concerned 
entirely with free fields, either as a model of the free electromagnetic field in an optical 
cavity, or as the zeroth-order building blocks in an interacting theory, e.g., the propagators 
as ingredients in higher-order Feynman diagrams. The development of the interacting theory 
and its application to CQED phenomena will be given elsewhere 0. A partial account of 
the present theory has been given in M. 



II. CLASSICAL FIELDS 

In this Section, we summarize the QNM expansion for classical fields. In this paper, we 
deal with scalar fields in 1 d only. 

For closed, linear systems, eigenfunction expansions, based on the eigenfunctions or NMs 
of their evolution operators, are a tool of vital importance in theoretical physics. However, 
open systems are not directly amenable to an NM analysis. Examples of open systems 
include optical cavities and finite regions of space near astrophysical objects, from which 
gravitational waves can escape In these systems, any initial state decays in time, so 

stationary NMs do not exist. As the simplest example, we shall be concerned with the scalar 
wave equation in one space dimension, 

pix)d^<P = dl^ (2.1) 

studied in a "cavity" < x < a, with the nodal boundary condition 

0(x=O, t) = 0(x=O, t) = (2.2) 

at one end but with the outgoing one 

0(a+,t) = -0'(a+,t) (2.3) 

at the other. The latter condition states that, just outside the cavity boundary, the field 
(j){x,t) is an outgoing wave 0(a; — t); the condition is specified at a"*" because, as we shall 
see below, one is often concerned with models in which there is a singularity in p{x) at 
X = a, leading to possible discontinuities in 0(a;) or 0'(x). The boundary condition ( p.3|) 
turns the cavity into a dissipative system that is leaky but not absorptive. The model ( p.l|) 
has been widely used as the scalar model of electromagnetism in an optical cavity . More 
physically, the 1-d nature is realized in Fabry- Perot cavities with lengths much smaller than 
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the lateral dimensions, and the scalar field model is rigorously applicable to the transverse 
electric sector. 

For the system ( p.l| )-( pl3| ), the eigensolutions, labeled by an index j, have the form ( p..!] ), 
with the QNMs or cavity resonances fj satisfying 

[dl + pix)co]]f,=0 (2.4) 

and the boundary conditions (|2.2|), ( |2.3| ) translating to 



/,(0) = 0, /;(a+)=zcu,/,(a+). (2.5) 

It is easily verified that Imuj < 0, so that the solution ( |1.1| ) is indeed decaying in time. 
Furthermore, the frequencies Uj, which we suppose to be ordered according to increasing 
real parts, are spaced by Auj Tc/a, approximately as for a conservative system of size a. 
With the possible exception of modes with Keujj = 0, the QNMs always occur in pairs with 
u_j = —ujj, and one can choose f_j = f*. 

The usual formalism concerning eigenfunction expansions relies on the hermiticity of the 
evolution operator, which only holds in the conservative case, and therefore breaks down for 
open systems. One possible resolution is to embed the cavity into a "universe" < x < A 
with a nodal condition at a; = A — > oo, and study its NMs — the modes of the universe. 
Namely, the system ( |2.1| )-( p73D is the restriction to x < a of the problem ( |2.1| ) on the half 
line < X < oo, if one sets 

p{x>a) = 1 (2.6) 

and with the extension of the initial conditions to the "outside" x > a obeying 
(j)'{x>a,t=0) = — 0(x>a, t=0). However, this has the obvious disadvantage of having to 
work with a continuum of states (spaced by Auj ~ vr/A 0) as opposed to the discrete 



set of eigenfunctions in the conservative case. Besides, the closed system of equations (|2.1|) - 
( p.3| ) shows that even in the presence of dissipation the time evolution of the cavity can be 
studied without explicit reference to the outside, which is the principal goal of the program 
of second quantization of the open system. 

Previous work (see [|10|-[T2[| and references therein) has established that, in spite of the 
lack of hermiticity in the conventional sense, an eigenfunction expansion for outgoing waves 
in classical open wave systems can be formulated in terms of the cavity degrees of freedom 
only, overcoming the disadvantages of the modes of the universe approach. The sufficient 
conditions for this QNM expansion are as follows. 



(a) The function p{x) has at least a step discontinuity at x = a. This demarcates a 
well-defined cavity region. 

(b) The function p{x) has no tail outside the cavity, i.e., p{x>a) = 1. This condition 
ensures that the outside does not refiect outgoing waves back into the cavity, enabling 
the complete elimination of the environment from the equations of motion. 

These conditions are satisfied for optical cavities bounded from extended vacuum by a sharp 
material interface. 
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The completeness of the QNMs can be pursued at two levels. First, one shows that the 
retarded Green's function of the system has the representation 



(2.7) 



for < x,y < a and t >0, where the /j's are normalized according to below. Thus, 



the dynamics is contained entirely in the QNMs, leading to a simple method of obtaining 
the retarded propagators and quantum correlation functions, as sketched in Section |T|. 

Second, realizing that the wave equation like any classical hamiltonian problem, 

requires both position and momentum to be specified as initial data, one introduces function 
pairs (f) = {(f), (p)^ with the conjugate momentum (f) = p(j), so that for eigenfunctions fj = 
{fj, —ipujfj)^. The space of all function pairs satisfying the boundary conditions ( p.2|) 
and (|2.3| ) will be denoted as F — the space of outgoing waves. 

Using these pairs, one can prove that the time evolution generated by ( |2.7| ) can be recast 



in the form [13 



(2.8) 



where the expansion coefficients are given by 

= 2^< 



(2.9) 



with aj{t) = aj{0)e and the bilinear scalar product 



{4>, X) =i{ dx 




[x)x{,x) + 0(x)x(x)l + 0(a+)x(a+) \ . 



(2.10) 



By simply letting t J,0 in ( p.8|) one arrives at a two-component expansion for an arbitrary (j) G 
F [|1^,|I^. This expansion makes the completeness of the QNMs manifest. The normalization 
used in 



to (|2.9| ) can be concisely expressed as 



/i' fj) - 2t^i 



(2.11) 



It is seen that ( |2.11|) in general is not real, underlining the difference between the product 
( p.lOj ) and a conventional one involving complex conjugation. The fact that (|2.11| ) is bilinear 
also serves to establish a phase convention for the wavef unctions. 
Upon introducing the two-component evolution operator 



H 



p{x) 



-1 
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(2.12) 



the cavity evolution ( |2.1| ) can be written as idt(t> = Hcf), in striking analogy with quantum 
mechanics. In this notation, the definition (|2.4|) of /-,• takes the form Hfj = ujjfj. The 
operator H can be shown to be symmetric with respect to the form ( |2.1U| ), i.e.. 



{<t>,Hx) = {x,Hcf>) 



(2.13) 



6 



for any (f),x ^ T. This analog of hermiticity holds even though the system is not conserva- 
tive. The symmetry of H yields the orthogonality relation 



(/,-,/a:)=0 foTUjJ^Uk (2.14) 

in an immediate transcription of the usual proof, leading to the uniqueness of the expansion. 
Incidentally, an expansion such as but involving the first component alone would not 
be unique. 

It is appropriate to contrast the two methods of approaching the classical theory, since 
they respectively underpin the two methods for dealing with the quantum case. The expan- 
sion of the Green's function is easy to derive, and is readily generalized to higher dimensions; 
however, in itself it does not lead to a unique expansion of the field, nor to concepts of or- 
thogonality. The two-component approach based on (|2.8| ) is more elegant, exhibits a deeper 
resemblance to conservative systems, and most importantly leads to a unique expansion 
in terms of orthogonal functions. The two-component expansion can in principle be gen- 
eralized to higher dimensions by treating each angular momentum sector / as a 1-d radial 



problem [16|, but the degree of complexity increases with /. Thus each method has its own 



merits; both will be pursued below, and the results compared. 

III. GREEN'S FUNCTION METHOD 

The quantum mechanics of the system is specified by the hamiltonian 

oo 



H = dx h(x 
Jo 




(3.1) 



together with the canonical equal-time commutation relation 

[<P{x),4>{y)] = i6{x-y) . (3.2) 

Time evolution is then generated by means of the Heisenberg equation A = i[H,A] for 
an arbitrary operator A. However, instead of the equations of motion for the quantum 
operators, in this Section we focus first on the retarded propagator 

G^(x,i/,t) = -z9{t)mx,t),^{y)]) , (3.3) 

in which is of course to be regarded as an operator, and (■ ■ ■) denotes the expectation 
value at a finite temperature T = 1//?; throughout we take h = = I. 

The central idea is that this propagator defined in terms of the quantum fields can be 
evaluated without explicitly introducing an expansion for the field operators (p, by simply 
noticing that G^{x,y,t) is exactly the same as the Green's function G of the classical wave 
equation [|1^, which has the expression (|2.7| ) in terms of QNMs. This relationship between 

and G follows from the commutation relation ( p.2|) . 

In terms of G^, it is straightforward to compute the equilibrium correlation function 

F(x,y,t) = (0(x,t)0(y)). (3.4) 
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We shall devote attention to F, because the physical quantities of interest in quantum field 
theory can often be formulated in terms of correlation functions, at either zero or finite 
temperatures. For example, the Casimir force is merely the vacuum expectation value of 
the electromagnetic stress tensor, which is an equal-time equal-space correlation function of 
two field operators. The spontaneous decay rate of an atom in an excited state is, in the 
golden rule approximation, related to the correlation function of two electric field operators. 
Since the correlation function is related to the retarded propagator W^, one gets 



F{x,y,u) 



iuj 



lmG^{x, y, uj) 

j2 fj(^)fj(y) 



The real-time correlator can be obtained from ( p.6|) by contour integration, yielding 



fjix)fjiy) 



-I3wj9(—t)—iujj \t\ 



+ E 



m=l 



G^{x, y, -ifXrn) - G^{x, y, i^r. 



(3.5) 
(3.6) 



(3.7) 



The first term in this formula is due to the QNM poles in F{x,y,Lj); the second term, 
which has no counterpart in G^{x,y,t), is caused by the Matsubara poles in F{x,y,u) at 
frequencies fim = 27imT. 

This very simple derivation has the advantage that it goes through in situations where 
the two-component formalism may be more complicated. 

In principle, physical quantities can be expressed in terms of F — bilinear quantities 
(such as the energy density) as linear combinations of F and its derivatives, and other 
quantities involving products of F's. For example, the energy density is 



{h{x)) = - \-p{x)df + d^dy] F{x,y,t) _ 

However, this quantity is divergent. Subtracting off the zero-point, we consider 

U{x,T) = {h{x)) - {h{x))T=o 

^ -p{x)dt + d^dy] Fs{x,y,t) 



(3.8) 



(3.9) 



in terms of the subtracted correlation function F^ = F — Fq, where Fq = lim^^oo-^- The 
limit /? — oo is best taken in ( p.6|) prior to Fourier inversion. We further make use of the 
expansion for G^ to get 



Fs{x,y,t)=Y.^^^^^C,{t) 



C,{t) 



3 



(3.10) 



/J oo fJ'TTli 



-^e^"^"*Ei(zcj,-t) - -^e-*'^^*Ei(-2cj,-t) , 



(3.11) 
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where this and subsequent formulas for Fs are written only for t > 0. Ei{z) is the exponential 
integral function [|T^] 



Ei(z) = / du , 

Jz u 



(3.12) 



in which the integration contour is defined not to pass through the origin and the negative 
real axis; on that semi-axis, the function is defined as the principal value. We further define 

Alternatively, for greater formal similarity to the conservative case, ( |3.11| ) can be rewrit- 
ten as 



1 ^ II p f^mt 



(3.13) 



(3.14) 



The prime on the sum in ( 3.13|) signifies that terms with HeuOj = are to be taken with 
weight |. 

The actual evaluation of U (x, T) needs some care in the j ^ oo part of the sums. These 
details, and the very similar problem for the calculation of the Casimir force, will be given 
elsewhere. 



IV. INCOMING WAVES 

The expansion of a classical field sketched in Section y is restricted to outgoing waves, 
i.e., to G r, satisfying ( p.3|) . In preparing the ground for the expansion of a quantum field, 
it is necessary to remove this restriction, for the simple reason that the zero-point quantum 
fluctuations will inevitably contain incoming waves as well. Moreover, one would wish that 
the ensuing theory should be applicable to situations where there is an incoming pump field. 

Thus, we study the wave equation ( |2.1| ) for the system together with the outside "bath" , 
i.e., on the half line x > 0, with p{x) satisfying (|2.6| ) and the boundary condition (|2.2| ). The 
initial conditions are now arbitrary and accordingly the outgoing boundary condition ( [^.3[ ) 
is abandoned, i.e., the restriction of cf) to the cavity need not lie in F. For the outside x > a 
(where p{x) = 1) the initial data are decomposed as 

cf){x>a, 0) = 0IN + 0OUT , (4.1) 

with satisfying the incoming wave condition = while (/>out = ~<^out- Fc>^ ^^e 
cavity subsystem this decomposition leads to the boundary condition 

0'(a+, t) + 0(a+, t) = 24>m{a + t) 

= h{t) , (4.2) 
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where the driving force b (see ( [4 .41 ) below for its name), being determined by the initial data, 
is supposedly a known function (at least in a statistical sense) which characterizes the waves 
incoming from the outside. Inside the cavity, the field is then expanded in terms of QNMs 
by (IH) with 




dx fj{x) (f){x) — ip{x)ujj(j){x) + fj{a'^)(j){a^) > . (4.3) 



That is, we retain the expansion formula and the inner product definition and notation 
even though (p ^ T. As a consequence, the sum in (p.8|) will in general not converge to 
at X = a, the point where the boundary condition is imposed. Nevertheless, the sum for 
the first component converges to everywhere, while the sum for the second component 
converges to everywhere except at x = a [|r^ . (This is most easily appreciated by noticing 
that upon changing at just one point, the resultant wavefunction can be made to lie in F.) 
This flaw on a set of measure zero does not lead to problems, however, for the projection 
formula ( [4.3|) renders the coefficients aj{t) well-defined irrespective of the convergence of the 
series ( p.8|) . 

The equation of motion for aj, which will survive quantization, will now be derived. By 
differentiating (|4.3|) with respect to time, and then integrating by parts, one obtains 



a,{t) + zuj,a,{t) = —f^{a+)b{t) . (4.4) 

In contrast to the case of pure outgoing waves, there is now an extra term on the right 
hand side: each QNM is driven by the "force" 6(t), and at the same time decays because 
of \m.ujj. The coupling to the "force" is determined by the surface value of the QNM 
wavefunction fj{a^). 



V. THE FIELD EXPANSION METHOD 



Another approach to second-quantization proceeds more explicitly by first promoting 
and to operators These fields may be regarded as operators for the entire "uni- 



verse" , which is a conservative system to which canonical quantization can be applied. The 
same projection formula (|4.3| ) as in the classical case now defines the a^'s as Hilbert space 
operators, obeying the equation of motion 



The crucial point is that the field commutation relation (|3.2|) and the projection formula 



now lead directly to commutators for these coefficients, viz., 

1 



dxdy [fj{x)fk{y)[(i){x),(j){y)] + fj{x)fk{y)[(j){x),^{y)] 
+ r''rfx/,(x)/fc(a+)[0(x),0(a+)]+ r''rfy/,(a+)/fc(y)[0(a+),0(|/)] 



JO 



(5.1) 
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In these equations, and are g- numbers, while fj, fk are c-cumbers. The two surface 
terms on the last line cancel as long as the delta function at the boundary of the integration 
interval is interpreted consistently. In the first line, the commutation relation (|3.2|) gives 
S{x — y) and cancels one integration. One is then left with 

AujjUJk 



[aj,ak] = '^ dxp{x)fj{x)fk{x) (5.2) 

4:UJjUJk •'0 

(5.3) 



where the second form follows from the first by means of the orthogonality relation ( p. 141 ), 
and will be useful later for comparison with results from Section 

The linear space structure for open systems involves projections based on the generalized 
inner product ( |2.10|) which is bilinear rather than linear in one vector and conjugate linear in 



the other; thus the expression in (|5.2| ) involves an integral over fj{x)fk{x) without complex 
conjugation. However, for the sake of a more transparent comparison with the conservative 
case, it is useful to re- write these expressions by changing j i— > —j and using a_j = aj, 
u^j = —u* and f-j = f* to give: 

[a],ak] = — ^ ^ I dxp{x)f*{x)fk{x) (5.4) 



z(a;; + a;fc)/;(a+)/fc(a+) 



(5.5) 



The result in the form ( |5.4| ) reveals the conservative limit most clearly; in this limit the 
integral would simply be 5|j|,|fci- 

Comparison of ( ^.41) and (|5.5| ) shows that 



^ 2 Imcj ' ^ ^ 



in the conservative limit. A more explicit proof is given in Appendix ^ 
The above commutators show that, if we define, for j > 



a] = \P^*a_^ , (5.7) 



then in the conservative limit these should reduce to the annihilation and creation operators, 
respectively Indeed, the QNM expansion ( |2.8| ) then takes the form 



j>0 



the standard normal-mode field expansion for a closed cavity [|T^. For finite damping, 
however, the operators have mixed creation and annihilation character. 
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In short, we have estabhshed an expansion of the quantum field (and its conjugate 
momentum 0) in terms of operators Oj and aj, and then obtained equations of motion 
and commutation relations for the latter. This, in principle, completes the program of 
second-quantization, and it remains to use these results to compute correlation functions 
and propagators, which we proceed to do in the following Sections. 

However, the deviation of the commutators ( |5.3|) and (|5.5|) from the canonical form 
prevents the construction of a Fock space, as is the case for quantum dissipative systems in 
general [0]. 



VI. CORRELATION FUNCTIONS 



The formalism derived in the last Section for expanding the quantum field (f) in terms of 
the operators at and aj will be applied to the calculation of equilibrium correlation functions, 
yielding discrete representations for the cavity correlator F in the presence of dissipation. 
Section [VI A| investigates the general case, and Section |VI B| compares the results with those 
obtained from the Green's function approach in Section |ITT| . Section |V1 C] evaluates and 
discusses the density of states. 



A. General case 

In equilibrium, the initial conditions for ( ^4.4| ) are irrelevant and the dynamics are com- 



pletely specified by the driving force b, i.e.. 



a^.(t) = f rft'e*-.(*'-*)6(t') . (6.1) 

The nonzero imaginary part of the Uj renders the integral rapidly converging, in contrast to 
the conservative case. Fourier transforming and taking expectation values then lead to 

{a,{u)a,) = ff^'-^f^^f^ (^(c.)6) . (6.2) 

Since b is fully specified by the mcoming waves from the free string a < x < oo, it does not 
"know" about the cavity x < a, so one can use the free infinite-string correlation function 
to calculate its spectral density from the definition (W^ as 



(b{u;)b) = -{d,-iu;Y{4>{x,u;)(f){y)),^Jy^^ 



2uj 



1 _ 



y=x 

(6.3) 



For a simple check, antisymmetrize ( |6.2|) in j and k and perform the inverse Fourier 
transform to reproduce ( ^.3| ) (for the expectation value of the commutator). Incidentally, by 
assuming other forms for {b{uj)b), the theory accommodates various incoming pump fields. 
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Given the two-point function for the driving force, it is straightforward to compute 
the two-point function for the response, namely the field-field correlation function inside the 
cavity. This now merely requires summation, that is, combination of ( p.8| ), ( |6.2| ) and ( |6.3| ) 
leads to 

F{x,y,u:) = —^Y. , ff^'^^f'^f^ . mUy) . (6.4) 



The above derivation leads to a clear physical interpretation of the pole structure of 
in the complex cu-plane: the Matsubara poles at = i^m = ^muiT {m G Z) arise from 
the thermal character of the incoming noise, while the QNM poles correspond to cavity 
resonances excited by this noise. 



B. Comparison of two forms for the correlation function 

It will be noticed that we have obtained two different QNM expansions for F, namely 
the double sum in ( |6.4|) and the single sum in (|3.6|) . We next prove their equivalence without 
invoking the QNM expansion of a quantum field. 

To do so, we rely on the identity p2 



G^{x, y, ij) - G^{x, y, -u) = — G^{x, a+, u) G^{y, a+, -u) (6.5) 

for x,y < a. This identity, proved in Appendix has no nontrivial counterpart in closed, 
conservative systems. For an interpretation, notice that G^{x,y,uj) — G^{x,y, —uj) oc 
lmG^{x,y,uj) vanishes in the conservative limit and hence is a measure of dissipation, 
which the right hand side states as taking place exclusively at the surface x = a"*". 

Given this identity, the equivalence of the two expressions for F follows simply by can- 
celing the Bose prefactors in ( p.6|) and ( |6.4|) and comparing the result with the Fourier 
transform of ( p.7| ). 

Although the two forms are equivalent, each has its own attractive properties. The di- 
agonal form ( |3.6| ) is simpler, while the non-diagonal form ( |6.4| ) is manifestly factorizable: 
F{x,y,uj) = A{uj)x{x,u!)x{y, —uj) [|^. Anticipating a similar structure for Feynman prop- 
agators, the non-diagonal form permits a quantum in one mode j to propagate to another 
mode k, while the diagonal form implies that the mode index is "conserved" . 

The expansion of correlations involving alone is not unique, on account of the doubling 
of QNMs {j and —j) compared to NMs [jll],|12|. As discussed in Section ||, it is more 



natural to consider (p = {(j), cj))^ , which leads to a unique expansion. Thus we define a tensor 
field-field correlator 

F(x,?/,u;) = {4){x,uj) ® 4){y)) 

^ ( {<P{x,u)<P{y)) {<P{x,u)^{y))\ 
\{4>{x,u)ct>{y)) {4>{x,u:)m)) 
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which can be expressed as 



f{x,y,uj) = J2(^jk{^)fj{x) ® fk{y) 

jk 

where ajk is evaluated to be (Appendix [C]) 

^fj{a''^)fkia+) 



(6.7) 



ajk[uj) 



2(1 — e f^'^)ujjUJk{ujj — uj){uJk + 00) 



(6.8) 



that is, the non-diagonal expansion (|6.4|) is the unique one which generalizes to the tensor F 
as in (|6.7[). 



C. Density of states 

Another important quantity is the density of states, which figures prominently in the 
heuristic argument of Purcell [^] and others [Q. The local density of states d{x,uj), given 
below only for real positive u, is related to the correlation function F by 

d{x,uj) = -—lmG^{x,x,uj) = -(1 -e-^'^)F(x,x,cj) , (6.9) 

TT TT 



which allows expression of this important quantity in terms of the QNMs. From (3^) one 
gets 

d{x, u;) = -y Im , , (6.10) 

while the non-diagonal expression ( |6.4| ) gives 



d{x,u) = —^ f —-fj[x)fk{x) (6.11 

= - E . ^^-^""If^^""^ , f,ix)Mx) . (6.12) 

The second form results from the first by use of the identity 

^^^•^^^^^•^^^=0, (6.13) 

which follows from ( |2.7|) by letting t | 0. 

Superficially, the diagonal form is simpler. However, if we take a single resonance ap- 
proximation, ( 6.101 ) yields, with one term j 



d{x, tu) ^ - Im /'■^ '-'^^ , , (6.14) 
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which is not positive definite. On the other hand, for the non-diagonal form, the appropriate 
approximation is to take one j and k = —j in ( 6.12|) , leading to 



\fM'-)fji^)\' 



^(^.^)- o,^..'T^.^f:..^,■.21 > (6-15) 



which is manifestly positive and moreover lorentzian. From this expression one finds, to 
leading order in |Ima;j| = 7, that 

duj / dx p{x)d{x,uj) ~ 1 , (6.16) 
5 Jo 



where the cj-integral is over one resonance. This statement is readily derived from ( |6.15|) by 



using ( |5.6| ) and the fact that Jdx p{x)\fj{x)\'^ ~ 1 for a narrow resonance. Recall that in the 
modes of the universe approach the unit weight of the resonances emerges simply as a 
numerical result, and is difficult to understand theoretically. Here the same result (in 1 d) is 
justified analytically, and moreover, one can in principle (a) estimate the corrections due to 
other resonances (note that there is no "background" apart from the QNM contributions), 
(b) calculate the corrections to higher order in 7, and (c) discuss the local density of states 
d{x,u}) rather than the integrated Jdxd{x,uj). Incidentally, this discussion shows that of 
the two equivalent forms (|6.11|) and (|6.12|) , the latter is the more appropriate, since it leads 



to a finite integral over u in the single-resonance approximation. 
One can derive another sum rule, 



d{x,uj') du' ^ — , (6-17) 
7vJp{x) 



for large tu. This second sum rule pj] states that the states are merely redistributed with- 
out changing their total number. However, this sum rule is not immediately useful when 
expressed in terms of the QNMs, and will not be further discussed here. 

VII. FEYNMAN PROPAGATOR 
A. Derivation of the Feynman propagator 

Another important correlation function is the Feynman propagator 

G^(x,y,t) = -2(T{0(x,t)0(i/)}), (7.1) 

in which T denotes time-ordering. Taking the Fourier transform of the definition ( [7. 1|) leads 
to a direct relation to the correlator (137 



G'^ix, y,co) = ~ + , ^ . I F{x, y, 00') . (7.2) 

We shall limit the discussion below to T = 0. Substitution of the right hand side of 
into (|73) yields & as fg^ 
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G^{x,y,u;) = -^Y^M^-}M^ U^)^^ + ^(-^)^| m Uy) . (7.3) 

2uJjUJk{UJj + LOk) [ iOj-LO Uk + UJ) 

The cavity Feynman propagator can also be expressed in diagonal form, by substituting 
( p.6|) into ( |7.2| ). Again taking T = 0, this leads to 

G^{x,y,u) = W^4rr^, (7.4) 
2 ^ ujj{\uj\ - ujj) 

2^M(|o;|-a;,) ' ^'''^ 

for real u. It is stressed that these forms as single sums exist even though {T{aj{t)ak}) 9^ 
Sj^±k in general. The form (|7.5| ) for has been derived from (|7.4| ) by means of the QNM 
identity ( |6.13|) . The second form with its divergence [26| at = is less convenient than 



the first. It has been included to show that caution is needed when speaking about "the 
contribution of one QNM". In fact the two summands are almost equal if \uj\ ^ Uj; such 
resonances are seen to be exclusively associated to terms with j > 0. 

All of these equivalent expressions (|7.3|), ( |7.4| ) and ( |7.5| ) can be written generally as 

G''ix,y,uj) = 5:/,(x)A,fcH/,(y) , (7.6) 

jk 

with different forms for Ajk- This has an obvious diagrammatic interpretation: the field at 
X (y) couples to the QNM j {k) with a vertex fj{x) {fk{y)), and the QNM propagates from 
mode j to mode k with an amplitude A^^. This may be compared with the more familiar 
case of an infinite conservative system, say 

G''{x,y,u) = f^e-'P^A{p,u;)e'^y . (7.7) 

It is seen that Jdp ■ ■ ■ is replaced by J2jk ' ' ' ■ An important goal of the present second- 
quantized theory is to study cavity-atom interactions 0, often referred to as CQED. The 
objective is to establish a set of "QNM Feynman rules", in which each line in a diagram is 
represented not by a continuous momentum, but by one discrete index (or a pair of them) 
— not only for computational convenience, but also because each term can be associated 
with a cavity resonance. Such a discrete representation is especially useful for microscopic 
cavities, where the resonances are widely spaced in frequency. The above results are crucial 
for establishing these Feynman rules. 

The possibility of alternate expressions for the propagator may recall a similar situation 



with gauge theories [2/], though the reasons are quite different. 



B. Decay rate and the resonance approximation 

While the use of the Feynman propagators in an interacting theory will be presented 
elsewhere |^ , it is neverthess profitable at this point to consider the very simple example of 
an atom coupled to the field at a fixed point x; in the dipole approximation, the decay rate 
is related to the equal-space propagator 
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(7.8) 



In particular, we shall be interested in the single- resonance approximation (RA) for D. 
The obvious choice is to take a single term of the sum in (|7.4|), i.e., 



D(uj) 



^ra' (t^) 



fj i^of 



2ujj{\uj\ - iUj) 



(7.9) 



The alternative is to start from ( |7.3| ), and retain only the (j, —j) + terms (only when 

k = —j does the factor uj+utk in the denominator of (|7.3|) get small close to the conservative 
limit, which is the only case in which a single resonance can dominate) to arrive at 



|/,(xo)/,(a+)p 



\LU\ 



UJj \UJ\ + iO* 



(7.10) 



Without loss of generality choosing j > 0, the second (non- resonant) term in ( |7.1(JD is of 
the same order as those already neglected, and hence for most purposes may be omitted. 
However, only the sum of the two terms in (|7.10| ) preserves the fundamental relation [jl7| 



D 



R, 



UJ 



(7.11) 



for real u, where {D^) is the retarded (advanced) propagator obtained from D{lj) by 
continuation from positive (negative) frequencies. As a consequence, it turns out that 
keeping both terms and using the ensuing cavity propagator to compute the self-energy of 
a two-level atom leads to a renormalization of the level splitting that is guaranteed to be 
real. Of course, for Dj-a' the equality ( [7.11| ) is always violated. 

Moreover, D^a' does not obey the equally fundamental inequality lmD{uj) < on the 
real axis [|T7| (see also ( |7.12| ) below), which D^a satisfies term by term. Violation of this 
inequality in general leads to a retarded atom propagator that has poles in the upper half 
cij-plane 0, signifying an unphysical instability. 

To be sure, in spite of these crucial differences between D^-a and Dra' their residues 
at |a;| = uoj agree in the conservative limit, in which the domination of a single QNM 
becomes rigourous. For a proof it suffices to note that |/j(a"'")p/2 jlmcujl — > 1 in this limit 
(Appendix 0). 

We have discussed the single-resonance approximation to both the density of states 
d{x,u) and to the equal space propagator D{lj). In fact, the arguments are equivalent, 
which can be appreciated physically from the fact that they both relate to the decay rate, 
and mathematically from the following identity for real positive u: 



d{x,uj) = lmD{uj) 

71 



(7.12) 



In several places we have remarked that the non-diagonal QNM representation has some 
nice properties, and is in fact the unique representation if the field and the conjugate 
momentum (p are considered together, for example in the tensor correlator ( |6.7|) . There are 
of course many ways to understand why the correlator is non-diagonal; one of the most direct 
is via ( |4.4| ), which shows that all the mode coefficients aj are driven by the same force b{t), 
so in general different coefficients will have phase coherence and hence a nonzero correlation. 
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Incidentally, this non-diagonal nature is not a quantum effect, since the property survives 
at high temperatures, e.g., /? — in (|6.4| ). However, all the propagators and correlation 
functions become diagonal in the conservative limit, as they should. In fact, applying ( p.l4| ) 
to ( |6.2|) in this limit readily yields 

{aj{u)ak) = ~ ^j)^j,-k , (7.13) 

in agreement with the creation-annihilation interpretation of the aj in this limit given above 
{ f).8\) . As a result both (|6.4| ) for F and (|7.3| ) for become diagonal in the conservative 



limit as well. 



VIII. EXAMPLE: THE DIELECTRIC ROD 



A useful check and example of the preceding is given by the "dielectric rod" model [10 : 



p{x) = n^9{a- x) +nl9{x - a) . (8.1) 



That is, we generalize the condition (|2.6| ) and allow p{x>a) to be an arbitrary constant n^. 
To be sure, this generalization is trivial in principle since a model with parameters (n, no, a) 
can be mapped onto one with parameters {n/no, l,noa) by the substitution x i— > UqX. Yet 
it is convenient in practice, since we can now deal with two different conservative limits (see 
the discussion below ( ^.4|) ): n/riQ (the "nodal" limit) and n/riQ oo (the "antinodal" 
limit), by letting no ^ oo and no — * 0, respectively, while keeping n, a, and hence Reujj 
(see ( |8^) below) fixed. 

The model ( ^.1| ) can be solved exactly for the QNM frequencies |]TT|, which read 



= J ij + - ^ acth(n/no), n > no ; 

I JTT — z ath(n/no), n < Uq ' 

1 , no + n , , 
= jTT - - In ^ . 8.3 

2 Uq — n 

Both no — i> oo and no — are indeed seen to be conservative limits {Imuj 0). They 
correspond to clamped and free ends, respectively, in the interpretation of the wave equation 
as the transverse vibrations of a string [ pU| , pU[ . On the other hand, for no — > n the QNM 
description breaks down as the dissipation tends to infinity. 
The QNM wavefunctions are given inside the cavity by 



/i(^) = \/i sininujx) ; (8.4) 



n^a 



their normalization is still given by ( |2.11| ), but our generalization p{x>a) = is readily 
shown to imply a corresponding modification of the surface term in the scalar product 
definition itself, viz., 



(0, X) = M / \<P{x)x{x) + ij){x)x{x) \ + no0(a+)x(a+) ) . (8.5) 
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In the "nodal" limit no —>■ oo, (|8.2|) and (|8.4|) show that fj{a) ~ rig , so that the 
factor fj{a)fk{a) in the surface term of the orthogonality relation (|2.14 ), (|8.5|) overcomes 
the explicit factor no, allowing the surface term to be neglected. This nodal limit has a 
counterpart in the "loaded string" model p{x) = 1 + M6{x — a), where M can be set to 
infinity ||31]|. On the other hand, in the antinodal limit tiq ^ it is the explicit uq which 
allows neglect of the surface term, fj{a) tending to a constant. Hence, the QNM expansion 
becomes a standard normal-mode expansion also if uq <^ n, which clarifies a detail left open 
in Section || and Ref. [|lT] . 

By means of a partial fraction expansion and the identity ( |6.13| ), the expression (^) for 
F can be rewritten as 



F{x,y,u) 



2cj(l 



UJ 



00 j + UJ 



(8.6) 



which is analytically more convenient even though the sum over j converges more slowly. 



Upon substitution of the dielectric rod QNMs (|8.2|) and (|8.4| ) into ( ^.6|) , the sum over j can 
be performed as the conventional Fourier series 



p-jiTTZ 

T.- — - 

111 — la 
3 •' 



1 ' 



implying 



^ jvr -ia~ - 1 



< ;z < 2 



and some rearrangement yields the correlation function as 



2nQ sin(nu;x) sm{nujy) 



UJ 1 



g-/3a;~) |^^2 siii^i^nuja) + 11? co??{nuja)] 
This result also follows from the modes of the universe approach in Appendix 0. 



(8.8) 



(8.9) 



As discussed in Section |T|, the subtracted correlation function Fq is directly related to 
the energy density, and is in fact the squared amplitude of the field strength. Figure shows 
Fs(x,a;,t) versus x aX. t = 0.1 for the dielectric rod model with a = 1, no = 1, n = 5, for 
different values of /5; this shows that the field amplitude is largest near the leaky end of the 
rod. Figure ^ shows -Fs(x, a;, t) versus t aX x = 0.3 (all other parameters the same as before). 
This diagram vividly illustrates the advantage of the QNM approach — although the result 
is in principle obtainable from the modes of the universe method, the clear oscillatory signal 
is best captured by expressing this quantity in terms of QNMs. 

Also, for the Feynman propagator the sum (|7.4|) can be performed if p = const, yielding 



G^{x,y,uj) 



sin(nc(jx) n cos[nco'(a — y)] — iriQ sm[n\uj\{a — y)] 



nuj 



ncosynuja) — mo sm[n\uj\a) 



.10) 



for X < y, while foTx>y the propagator is obtained via G^{x,y,u) = G^{y,x,uj). 

Notice that the final expressions ( p.9|) and ( 8.10| ) tend to a finite limit if ng — >■ n even 
though the individual terms in ( ^.6] ) and ( [7. 4] ) do not. In this semi- infinite string limit the 
very notions of cavity and environment lose their meaning, and indeed the right hand sides 
of ( ^.9|) and ( p.lO|) are seen to become independent of a. 
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IX. FINAL REMARKS 



To summarize, we have developed the second-quantized version of the field theory using 
the QNM basis. Various physical quantities are then written as sums over QNM contribu- 
tions — either as diagonal sums over a single index j, or as non-diagonal sums over a pair 
of indices jk. The resonance approximation is studied, leading to a proof of the unit weight 
of narrow resonances in the density of states, or equivalently the enhancement rate for the 
decay of excited states as embodied in the behavior of the equal space propagator D^u). 

As has been mentioned already in Section |I|, an important extension of the present 
work is to include matter in the hamiltonian enabling the application of QNMs to 

quantum optics. This will be the subject of the sequel paper Ref. [If. Other generalizations 
include the study of vector fields, and of open systems in three space dimensions. Further, a 
development parallelling the present one could be carried out for the Klein-Gordon equation 
instead of the wave equation ([2.11). Since the two evolution equations are directly related by 
a transformation of the spatial variable [|1T| however, this has not been taken up here. 

Instead of generalizing the physical system one can also relax the assumption of global 
equilibrium made in Section [V^. It is recalled here that the formalism of Sections ^ and [V| 
— and in particular the driving force b of (|4.2| ) — is well-defined for any initial state of the 
fields; taking a coherent state for the latter instead of a thermal one enables the study of a 
pumped cavity. 

On the theoretical side, it would be interesting to provide a path-integral formula- 
tion of QNM quantization. This can supposedly be done on two levels. The first, semi- 
phenomenological one is to write down an effective action generating dynamics equivalent 
to ([4.4] ). The second, more fundamental one is to start with the action for the whole universe 
for our model ( p.l|) , integrate out the degrees of freedom of the outside, and use a QNM 
basis for the ensuing dynamics of the cavity. 

For such future developments, this article can hopefully serve as a starting point and 
reference. In conclusion, we have shown that the QNM expansion is as powerful for open 
second-quantized systems as it is for their classical counterparts. 
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APPENDIX A: RELATION BETWEEN SURFACE TERM AND IMAGINARY 

PART OF THE FREQUENCY 

In this Appendix we give an alternative proof of the identity ( |5.(j| ) in the conserva- 
tive limit. For this purpose, generalize to complex classical fields and define the en- 
ergy density h{x) = \dx(p\'^/2 + |0p/2p(a;), so that h{x) = —dxj{x) with the current 
j{x) = — He (J){x) dx(f){x)^ /p{^)- Define the cavity energy E = jQ^dx h{x), then for a 
field (f){x,t) = /j(x)e~*'^^* at t = one has 
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E = ]-\f,{a+)\' + {Reu;j)' f^x p{x)\f,{x)\\ (Al) 
I Jo 

where 7 = llmcUjl. In the conservative hmit the first term vanishes, while the integral in 
the second term tends to unity so that E —>■ \ujj\'^. Combination with —E = 2'jE = j{a'^) = 
|u;jp|/j(a"'")p shows that \fj{a'^)\'^/2'j — > 1 in this limit, which proves our assertion. 



APPENDIX B: IDENTITY FOR RETARDED PROPAGATOR 



In this Appendix we derive the Green's function identity ( |6.5| ). To this end, define /(x, u) 
{g{y,uj)) as the solution of the homogeneous wave equation ( p.4| ) (upon the substitution 
ujj ^ uj) satisfying the first (second) of the boundary conditions ( p^.S] ) P^ . This allows one 
to write 

G-(.<y,.) = li^^l^^ , (Bl) 

where one can choose f{x,uj) = f{x, —u) = f*{x,uj) and g{y,u!) = g*{y, —uj), and where W 
is the position- independent wronskian of the functions / and g |Tn]. Then one has 



G'^ix, y, u) - G'^ix, y, -uj) g{y, uj)W*{uj) - g*{y, u;)Wiu;) 

(B2) 



G^x, a+, uj) G^y, a+, -uj) \g{a+, uj)\^f{y, u) 

9{y, ^)9'*{y, ^) - g'iy, ^)9*{y, ^) 



\9ia+,u;)\' 

The numerator of this last expression is itself another wronskian, and hence can be evaluated 
at ?/ = a"*" to yield 

G''ix,y,uj)-G''{x,y,-uj) ^ u;) 
G^x, a+, u) G^{y, a+, -uj) ' g{a+, u) 

_ 2uj 
I 

completing the proof of (|6.5|). 



(B3) 



APPENDIX C: EXPANSION OF TENSOR CORRELATOR 

The coefficients cijk in (|6.7|) are given by the projection 

. «!M^ (CD 

in terms of the bilinear form on the product space, which reads 
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((P>Q)) = - dxdy {Pii{x,y)Q22{x,y) + Pi2{x,y)Q2i{x,y) 
Jo 

+ P2i{,x,y)Qi2{x,y) + P22{x,y)Qii{x,y)} 
- dx \^Pii{x,a^)Q2i{x,a'^) + P2i{x,a^)Qii{x,a^)^ 

~ lo {^ii(a^y)<5i2(a^y) + ^i2(a^y)Qii(a^,y)} 
-Pn(a+a+)gn(a+ a+) . (C2) 

Substitution of ( |6.6| ) for F(a;) into (Pl|) yields 

djk{uj) = — l{uj +uj){ujk-uj) I dxdy p{x)p{y)fj{x)fk{y)F{x,y,uj) 



+i{ujj + uo)fk{a^) / dx p{x) fj{x)F{x, a'^,uj) 







+i{uJk - uj)f){a+) / dyp{y)fk{y)F{a\y,uj) 

•J 

-/,(a+)/fc(a+)F(a+a+u;)j. (C3) 



Inserting any QNM expansion F{x,y,uj) = J2imbim{^^)fi{x)fm{y) such as ( PTB| ) or ( p^ ) and 
invoking the relations ( |2.11|) , (|2.14|) and (|6.13|) , the expression (|C3|) can be evaluated as in 



which is what we set out to show. 

APPENDIX D: MODES OF THE UNIVERSE APPROACH TO THE 

CORRELATION FUNCTION 

It is instructive to rederive the correlator F using the modes of the universe (MU). The 



MU expansion of the fields reads |^| 



i.e., it is of the same form as (|5.8|) but in ( |D1|) the sum runs over the MU frequencies 
i>i = Itt/A (to leading order in a/ A <^ 1), the ui and ipi being MU annihilation operators 
and wavef unctions, respectively. Insertion of ( pl| ) into (|3^ ) yields 



Fix, y,uj) = }_^ {{ul (uj) + Ui{u)}{ul + u^}) 

= 1^ 7> W i^)ui + ui{uj)ul) 

I ^^l 

= ^Mx, \u\)ij{y, \uj\){e{-uj)N{-u) + e{u)[N{u) + 1]} 

\UJ\ 

^ AV^(a:, \uj\)i){y, \uj\) 
uj{l - e-P^) 



:d2) 
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To arrive at the third hne, we used u\^\uj) = 27rS{uj W i'i)u\^^ and defined the boson 
occupation number N{uj) = [exp{Puj) — 1]^^. Comparison of ( |6.4| ) and ( |D2| ) elucidates 
why the former factorizes with respect to x and y: this is seen to be a consequence of the 
nondegeneracy of the MU spectrum of the semi- infinite string, as opposed to, e.g., a free 
string or one with periodic boundary conditions. Evaluation of for the "dielectric rod" 



of Section [VII1| at once shows that ([D^) indeed reduces to ( ^.9|) in this case, providing a 
further comparison between the various techniques of this paper. 

Finally, in terms of the MU the local density of states is defined as 



d{x,uj) = ^ z/i)p(5(z/i -uj) 
I 

A ,2 



TT 



(D3) 



and comparison with ( P2|) at once reproduces (|6] 
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FIGURES 



FIG. 1. Equal-space correlation function within the dielectric rod as a function of x at t 
and different inverse temperatures (3. The refractive indices are no = 1, n = 5; the width a - 

FIG. 2. Equal-space correlation function within the dielectric rod as a function of t at x 
and different inverse temperatures /3. The refractive indices are no = 1, n = 5; the width a - 
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